#! /usr/bin/env python

import pyx
import pyfeyn
from pyfeyn.user import *

fd = FeynDiagram()
processOptions()

in1 = Point(-6,  2)
in2 = Point(-6, -2)
in11 = Point(-2,  2)
in22 = Point(-2, -2)

out22 = Point(4, 0)

vtx = Vertex(1, 0, mark=CIRCLE)

Gluon(in1,in11).addLabel(r"\Pgg")
Gluon(in2,in22).addLabel(r"\Pgg")
Fermion(in22,in11).addArrow().addLabel(r"\Pqt/\Pqb")
Fermion(in11,vtx).addArrow().addLabel(r"\Pqt/\Pqb")
Fermion(vtx,in22).addArrow().addLabel(r"\Pqt/\Pqb")
Higgs(vtx,out22).addLabel(r"\PHiggs")


fd.draw("./GluonFusion_Higgs_production_Feynman.pdf")

